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ABSTRACT 


A normal star is a mass of compressible gas of low Prandtl number, 
held together by its own gravitational force. It generates energy 
principally in its central regions, either by very slow gravitational 
contraction or, more commonly, by thermonuclear reactions; this energy 
is subsequently transported to the surface down a temperature gradient 
by radiative diffusion or convection. 

If the star has no angular momentum the equations governing its 
motion admit spherically symmetric stationary solutions. But if the 
star is rotating it is not possible in general for it to exist in a state 
with no motion other than its rotation; meridional currents are set up 
which advect angular momentum through the star. The details of this 
process are not well understood. Stability calculations seem to imply 
that no stable steady solutions exist. 

The distribution of angular momentum in a star depends on a 
balance between transport by the large scale meridional motions and by 
the Reynolds stresses of turbulence. Neither theory nor observation- 
is yet sufficiently refined to deduce what that distribution is even 
in the sun* 


INTRODUCTION 

The structure and evolution of the majority of stars for most 
of their life seem to be well represented by theoretical models that 
ignore the fact that they rotate. This is because the ratio of the 
centrifugal to gravitational forces, or rotational Froude number, is 
low, and for many purposes the dynamical effects of rotation produce 
only a small perturbation from a nonrotating configuration. Nevertheless 
rotation presents some very interesting fluid dynamical problems which 
are the subject of this review. 

The as tr ©physical consequences of rotation in stars are not 
without interest. Rotation affects the spectra of the emitted radiation 
from which conditions in stellar atmospheres are deduced, and the possibility 
of rotationally induced mixing, which is discussed below, has important 
implications concerning the nuclear transmutations that occur in stellar 
cores. In addition, there are stars which appear to be rotating so 
rapidly that the dynamical effects of the centrifugal force have a 
significant influence on their evolution. These matters are not discussed 
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here; instead the reader is referred to the reviews by Strittmatter [1] 
and Ericke and Kippenham [2] , and the IAU colloquium proceedings edited 
by Slettebak [3J . 

Stars condense from the interstellar medium and at some point 
attain a state of hydrostatic equilibrium* They are opaque and hot 
in the centre, the stored heat maintaining a pressure gradient which 
balances gravity . Their structure is such that the balance 

2 T + 1/ - 0 (1*1) 


is satisfied, where V is the gravitational potential energy of the 
star and T the total kinetic energy, both macroscopic and microscopic 
[43. This is the virial theorem of Poincard and Eddington, and assumes 
that no body force other than gravity is exerted on the star. Normal 
stellar material is gaseous and can be assumed to satisfy the usual 
equations of fluid dynamics* The perfect gas law 



( 1 . 2 ) 


is approximately satisfied, where p, p and T are pressure, density and 
temperature, R is the gas constant and u is the mean molecular weight 
of the gas (mean mass per particle measured in units of the hydrogen 
atomic mass) which depends upon the chemical composition and the state 
of ionization* Radiation pressure is usually comparatively small 
and will be ignored in this discussion* 


If dl the kinetic energy of the stellar material is in microscopic 
particle motions (so there is no fluid motion) equation (1.1) may be 
written 


3(y - i) a + v m o. 


(1.3) 


where (i is the total internal energy of the gas and y an appropriate 
mean ratio. of principal specific heats. The total energy of the 
star is 

E * a + v * - (3y - 4 } a. (1.4) 


It can be shown that for the star to be dynamically stable y must 
exceed 4/3. Thus as energy is lost by radiation from the stellar 
surface, E decreases and U increases. The star contracts slowly, 
liberating gravitational energy which is partly radiated away and 
partly stored in the compressed heated gas. This process, first 
discussed by Kelvin and Helmholtz, is controlled by the rate radiation 
can escape by diffusion from the star. It continues until the centre 
becomes so ho t that_nuclear reactions begin to take place. Contraction 
halts; the star has reached the main sequence, the total energy output, 
or luminosity, at the surface being supplied entirely by the thermo- 
nuclear reactions. Its gross structure now depends primarily on 
just its mass. On the main sequence stars convert hydrogen, their 
principal constituent, into helium. This is the longest phase of 
their evolution, and consequently most of the visible stars in the sky 
are main sequence stars. 
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Once hydrogen has been exhausted from the reacting core further 
gravitational contraction and additional nuclear reactions take place, 
the details depending on the star’s mass and composition, and eventually 
the star contracts into its final state : a cooling white dwarf supported 
against gravity by a sea of degenerate electrons, a neutron star of 
degenerate nuclear matter, or perhaps a black hole* 

This review does not attempt to discuss the properties of a 
star in all stages of its evolution, but concentrates on the main 
sequence* Magnetic fields, which are undoubtedly important in some 
stars, are hardly considered; it seems prudent to try to understand 
first what is probably the simpler dynamics of a nonmagnetic star. 
Observational techniques are not discussed at all. 

If one plots the mean surface rotation rates of main sequence 
stars as a function of their mass M the most striking feature is an 
abrupt rise, in the vicinity of 1| solar masses, from approximately^ 
zero to a value corresponding to a Froude number which is an appreciable 
fraction of unity [53 . Observations of the distribution of angular 
velocity across the surface have been made only for the sun, which is 
a typical slowly rotating star. The angular velocity 0 varies by about 

3 2 

20% from pole to equator, the rotational Froude number R # H /Qf*being about 

2 x 10 ^ £63 , where M 0 and R« are the solar mass and radius and G the 
gravitational constant. These are the principal main sequence observations 

* 

a theory of stellar rotation must explain. 

Aside from the as trophy sical and fluid dynamical aspects of the 
solar rotation interest has recently been aroused in connection with 
gravitation theory. One of the major puzzles of the last century was 
the discrepancy between the measured precession of the perihelion of 
the orbit of Mercury and the predictions of Newtonian gravitation 
theory [73 • One of the proposed explanations was that the sun was 
oblate, but the observed oblateness [8,9] turned out to be too small 
to account for the precession [103. With the advent of General 
Relativity, however, it was generally agreed that the problem had 
been resolved, because Einstein’s theory predicted the observed result 
[11] . But more recently Dicke [123 and Roxburgh £133 proposed that 
the apparent agreement was fortuitous, and that an alternative theory 
of gravitation [14] based on the scalar-tensor theory of Jordan £15] , 
which predicts a lower precession rate, is correct, the residual 
precession of Mercury’s orbit arising from an assumed oblateness of the 
sun’s figure. Subsequently the shape of the sun’s image was remeasured 
by Dicke and Goldenberg [16, 17]. An oblateness of about the 
required value was reported, if it is assumed that the oblateness 
of the image implies a similar oblateness in the shape of the gravitational 
equipotentials . That this assumption is justified has been argued 
by Dicke but has not been generally accepted [17,18]. This important 
issue will be discussed in the concluding section of this review. 


* 

Perhaps the most important question of the subject is why stars exist 
at all, for had they condensed conserving angular momentum to their present 
size from interstellar material with initial rotation rates typical of 
the galactic rotation they would not be gravitationally bound. This 
is a problem in star formation. 
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In the next section possible steady solutions of the equations 
governing rotating stars are discussed. Their stability is considered 
in the subsequent section, and the final section is addressed to how 
angular momentum is distributed thoughout a star as it evolves. But 
before that a few remarks about the magnitudes of quantities pertaining 
to stars are perhaps not out of place. 

The sun, which may be considered representative, has a mass 
33 10 

M # « 2 x 10 g, a radius R « 7 x 10 cm and a luminosity 
33 -1 

L # » 4 x 10 erg sec , At a distance of about 0.5 R # from the 
centre its density is about the same as its mean density which is 
approximately the density of water, and its temperature is about 
3 x 10 6 °K; the central density and temperature are about 100 g cm ^ 
and 1.5 x 10^ °K. At the visible surface, or photosphere, its density 

and temperature are about 3 x 10 7 g cm 3 and 5800 °K. The outer 
20% by radius is in a state of turbulent convection and the rest is 
generally believed to be convectively stable. Qualitatively this is 
a property of all stars at the low mass end of the main sequence 
Glower main sequence}; the higher mass rapidly rotating stars of the 
upper main sequence have only thin weak convective regions near the 
surface, and they also have convective central cores. 

Energy is transported through the convectively stable regions of 
a star principally by radiative diffusion. Except near the surface, the 
energy flux F satisfies 

F - - KVT , (1.5) 

3 

where K * 4acT /3 X p is the radiative conductivity; a is the radiation 
constant, x the opacity of the gas and c the velocity of light. At a 

median point in the sun its value is about 10^ erg cm~* sec~* °K~^. 

The kinematic -shear viscosity arising from photon momentum transport 
4 2 2 *1 

is « 4aT /15 xP c s 1 cm sec at a radius r of 0.5 R # and increases 

with r nearly everywhere. Viscosity due to microscopic particle ion 

notions v. is of order m p * (kT) 5/2 / (e*p InA) where A - (o p /p) (kT/e 2 ) 3 , 

m and e are the proton mass and charge and k is the Boltzmann’s constant; 
p . . 2-1 

this is about 10 cm sec at r » 0.5 R* and increases slowly with r. 

From these values can be calculated a timescale characteristic 
of large scale dynamical motions : the time for free fall under gravity 
from surface to centre. In view of the balance (1.1) this is of the 
same order as the sound travel time t and is about 1 hr. The thermal 
diffusion time across a distance R 0 is P equal to the Kelvin-Helmholtz 

gravitational contraction time t_^ = 10^ yr and is short compared with 

the nuclear time ^ 10 yr to convert say 10% of all the hydrogen 

to helium at the current luminosity. The age of the sun t # is 5 x 10^ yr. 
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The Reynolds number associated with the smallest velocity that can be 

2 

of global significance in this time, R*/t* , is Rq /(v fi +vpT 0 * 3000* In 
other words, the characteristic viscous diffusion time rs much larger 
than the age of the sun* For flow velocities comparable with the 

rotation speed (about 2 km sec ^ at the equator) the Reynolds number 

is of order lO 1 ^ throughout most of the star. Under most circumstances, 
therefore, viscous forces can be ignored. Only in the extreme outer 
regions and in boundary layers might this not be so . The Prandtl number 
(Vr+v^) pCp/K, where is the specific heat at constant pressure , is about 

* 10 ^ at a median point. 


STEADY ROTATING CONFIGURATIONS 

As with most physical systems it is expedient to seek the steady 
configurations first. They are comparatively simple to analyse and, 
if stable, are states towards which one might expect stars to evolve. 

In most of the published work steady axisymmetric models are 
constructed in which the only motion is a rotation about the z axis. 

The momentum equation then reduces 

- Vp * p7$ + pQ 2 m » 0, (2.1) 

where is is the cylindrical radius vector from the z axis, 0 is the 
angular velocity and $ is the gravitational potential which satisfies 

V 2 $ - -4irGp. (2.2) 

In equation (2.1) viscous forces have been ignored and it has been 
assumed that no magnetic field is present. 

The simplest models to analyse are of white dwarfs. In these 
stars the electrons are degenerate and a bary tropic equation of state 
is a good approximation. Because the fluid is assumed inviscid the 
angular velocity distribution is not determined completed by the 
equations of motion. However, it follows immediately from equation 
(2.1) that the centrifugal force is derivable from a potential V and 
that 0 varies only with distance IS from the axis of rotation. Accordingly 
equation (2*1) reduces to 


Vp - pV(| * V) 2 pW, (2.3) 

from which it follows that p and p are constant on surfaces of Hf , 
or level surfaces. There remains considerable freedom in the choice 
of 3(tff) • It can be set by specifying, for example, the angular momentum 
per unit mass h(m) , where m is the mass enclosed in a cylinder tsj * constant. 
Then any choice of h(m) leads to a solution. It was pointed out by 
Ostriker and Hark (19] that G(tff) cannot be specified arbitrarily. 

An important consequence of rotation concerns the possible masses 
of white dwarfs. The equation of state of the degenerate white dwarf 
material is such that in the special case of no rotation there is a 
maximum mass H, which can be supported against the self gravity of the 
star. This was first demonstrated by Chandrasekhar [20] who found 
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M *1.44 H 9 . The estimate has been revised with the use of a more 

refined equation of . state [e.g. 211, but the principle is unchanged and 
for a long time it was commonly believed that M c was an upper limit 

to the masses of real white dwarfs- However, it was emphasised by Mestel 
[22], for example, that this limit applies only to nonrotating spherically 
symmetric configurations, and Ostriker, Bodenheimer and Lynden— Bell [233 
suggested that rotating white dwarfs with masses much greater than M c 

exist- Because gravitational forces scale with their characteristic 

length L as t and centrifugal forces as L for a star whose mass and 
angular momentum are conserved, a rotating star more massive than M c 

would be prevented from collapsing indefinitely by the centrifugal forces. 
Indeed Ostriker et al- pointed out that were it not for the changes in 
physics that occur at extremely high densities a strict upper bound to 
the possible masses of rotating white dwarfs would be removed entirely, 
whatever the angular momentum. 

More interesting from a fluid dynamical point of view are the 
nondegenerate main sequence stars for which temperature is present in 
the equation of state. In these stars energy is transported down the 
temperature gradient from centre to surface- How configurations in 
which the only motion is rotation with an associated conservative 
centrifugal force field do not occur, and much of the wisdom gained 
from the simple barytropic models is not very useful- In general the 
energy transport equation does not permit T to be constant on level 
surfaces. Consequently, in view of equation (1-2) , p and p cannot both 
be constant on level surfaces, which implies that the momentum equation 
(2-3) cannot be accurately satisfied- The unbalanced pressure gradients 
drive a circulation u in meridional planes which advects both heat and 
angular momentum through the star [24-26] • 

Except in the surface layers of the star the motions induced 
are so slow that inertial forces are negligible compared with the other 
terms in the momentum equation so that equation (2-3) is very nearly 
satisfied- For a given rotation field the circulation velocity can 
thus be estimated from the energy equation 


pC p u.TT - 6 u.7p » pe - div F (2-4) 

simply as that required to transport the requisite amount of heat to 
maintain p, p and T (very nearly) constant on level surfaces. In 

equation (2.4) — g® and e is the thermonuclear energy generation 

P 

rate per unit mass- Thus equation (2.5) may be rewritten approximately : 


PTC 

in-ail 

P 


(7 - 7 &d ) u.7$ 


£ - i div F 
P 


(2.5) 
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where 7 * and V ad v frlnp / * tiie derivative being taken at constant 

specific entropy s, and 

div F - - f'(f)(?$| 2 + f($)[i ~<tff 2 n 2 ) - 4ttGo] , (2.6) 

jm 

where f ({) « K — and a prime denotes differentiation with respect 

to the argument. It is immediately clear from (2.5) and (2.6) that 

u cannot be zero if ft is constant, for since all the terms in (2.6) 

except 7$ are then constant on level surfaces it would follow that 

f f ($) * 0 and e * 2 (ft /p - 2irG)f j in general this cannot be satisfied 
because s depends on nuclear physics, and hence on p and T but not 
on ft, and in particular is essentially zero in all but a central core 
of the star [27} • There is no general proof that u - 0 when ft is 

a function of tr alone, though it has been shown to be the case if 

~ -g™ (tr 2 ft 2 ) - 4irGp vanishes in a region in which e m 0 [281 . Detailed 

calculations of the meridional flow were first performed by Sweet [29] 
who assumed slow uniform rotation. 

The magnitude of thie meridional Eddington-Sweet circulation 
velocity u can be estimated from equations (2.5) and (2.6) to be [30] 


u 




(2.7) 


where M, L and p are the mass, luminosity and mean density of the star, 

2 

g is gravitational acceleration and X « Eft /g is the Froude number. This 
estimate has been confirmed by detailed computations for certain 
specified nonuniform rotation laws with A. small [31] by expanding in 
powers of X about the nonrotating configuration; the case of uniform 
rotation is singular but appears to yield similar velocities [29, 32] • 


So far it has been implicitly assumed that the stellar fluid 
is chemically homogeneous. Nuclear transmutations in the cores of 
main sequence stars generate gravitationally stable gradients in the 
mean molecular weight p. According to Mestel [33] these prevent the 
Eddington-Sweet circulation from penetrating the core by deflecting 
them in a thin viscous boundary layer. 


It has been of some concern that the Eddington-Sweet circulation 
velocity, in the present approximation, diverges at the surface of the 
star where the density is almost zero. This arises because the heat 
capacity of the gas per unit volume is too low to transport at low 
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velocities the flux of heat required to satisfy equations (2.5) and (2.6). 
Smith [34] pointed out that excessive demands were being made by the 
diffusion approximation (1.5) near the stellar surface, where it isn’t 
valid, and showed that the situation was alleviated somewhat by using 
a more careful treatment of radiative transfer for the outer layers [see 
also ref. 353 * The velocities were still so great near the surface 
that the assumptions of the procedure, such as the neglect of the inertia 
terms, are hardly plausible,' and Smith concluded that in any case the 
flow would be liable to shear instabilities. 

Smith's analysis, in common with earlier work, assumed an inexorable 
uniform rotation, a weak poloidal magnetic field having been invoked 
to preserve it. In the absence of such a field, advection of angular 
momentum by meridional circulation is inevitable, and a state of nonuniform 
rotation would rapidly be set up in the surface layers. Osaki [353 
argued, without proof, that in that event readjustment of the angular 
velocity distribution would be such as to reduce the meridional circulation 
velocity in the stellar atmosphere to zero. 


Another place where the circulation is unable to advect heat 
efficiently is at the edge of a convection zone where V this time 

because the slowly rising or falling fluid, which maintains- pressure 
balance with its surroundings, cannot modify the temperature difference 
between it and its environment. Once again equation (2.5) demands an 
infinite velocity unless div F » pe when V * V^, which is not generally 

the case for an arbitrarily specified distribution of angular momentum. 
Osaki [373 pointed out once again that this implies that a steady solution 
does not exist for the assumed rotation law, and that appropriate 
readjustments must take place. 


From (2.7) can be estimated the timescale for stellar material 
to traverse a distance comparable with the radius R of the star : 


BS 


- CK 


KH 


- 2 * 


( 2 . 8 ) 


For the sun * 100 if it assumed to rotate with about its surface 

angular velocity throughout, which suggests that little advection of 
angular momentum by Eddington-Sweet currents is taking place. This 
conclusion must be treated with some caution, however, because the 
estimate (2.8) is very approximate and there is evidence [381 that it 
overestimates the circulation time somewhat. For the rapidly rotating 
. -2 

upper mam sequence stars t_* 10 t and considerable readjustment 

£b n 

of the rotation field must have occurred. 


The advection of angular momentum by meridional currents has 
not been discussed in great detail, though- Sakurai [383 has considered 
the early stages of evolution of a simplified solar model from a state 
of rigid rotation. A steady state with circulation must have 
2 

m 0 constant on streamlines, which implies that 0 is generally greater 
near the rotation axis. However, Mestel [303 has suggested that a 
star that conserves its angular momentum will redistribute it in such 
a way as to choke off completely the meridional flow. Accordingly one 
might seek (nonconservative) rotation laws which drive no Eddington-Sweet 
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circulation. Schwarzschild [39] , Roxburgh [40,41] and Clement [42] 
have attempted this, but their models are not entirely consistent. 

However even_ if_conaistent^4ao4e-ls— obtained , either with or without 

circulation, it is not at all obvious that a star could actually evolve 
to such a state. Furthermore, as is discussed in the next section, 
it seems likely that all steady solutions are unstable. 

This discussion has concentrated on the rotation of regions that 
are stable to convection. In convection zones too meridional currents 
are likely to arise except for very special distributions of angular 
momentum [43] > but here the dynamics is complicated by the influence 
of the rotation on the convection. It is common practice to assume 
the Reynolds stresses act to force rigid rotation [39*42] , though other 
assumptions have been made when discussing the convective envelope of 
the sun [see final section] . 

There has been no detailed study of nonaxi symmetric configurations 
of the type discussed by Lebovitz [44] using an equation of state that 
is realistic for stars. 


STABILITY 

The discussion will be limited to instabilities that arise directly 
from the rotation of the star. Dr Lebovitz has already explained what 
is known about fission [44] and I shall dwell on it no further. Of 
the modes of instability that remain the most widely studied have been 
the axisymmetric ones on a lepgth scale small* compared with the characteristic 
scale of variation of the basic ’equilibrium* state, and which derive 
from the Rayleigh instability.- 

Rayleigh [45] argued, by analogy with the instability of a 
density stratified fluid under gravity and by energy considerations, 
that a uniform incompressible inviscid fluid, steadily rotating between 
two coaxial cylinders, is stable to axisymmetric perturbations if 


4 C*V> o 


(3.1) 


everywhere in the fluid, and is unstable if (3.1) is violated anywhere. 

For an inhomogeneous incompressible fluid the criterion is 

4 (pA 2 ) > 0. (3.2) 

It is possible to perform a proper perturbation analysis to 
establish that criterion (3.2) is both necessary and sufficient for 
stability to axisymmetric perturbations [46] . Generalization to the 
more complicated flows encountered in stellar models is difficult, however, 
and one must resort to other arguments. A powerful method which is 
relatively simply to handle is to use the principle of virtual work ; 
an infinitesimal virtual adiabatic displacement £ of the fluid is 
imagined to occur and the energy difference 6E ~ between the final 
and initial states is computed. If <5E > 0 for all possible £ 
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there is no energy available in the fluid to drive any disturbance, 
and the basic state is stable. On the other hand, if displacements 
£ exists for which 6E < 0 instability cannot be deduced (unless 


5E < 0 for- all- §)- unless—it-a an be d cmo nstr ated~Xhat 6E is negative 
for a realizable disturbance which satisfies the equations of motion. 
However, even though instability cannot be strictly proved, in most of the 
special cases that have been investigated instability has been found. 

The method has been used by H0iland [47] who used d 1 Alembert 1 s principle 
to estimate the frequencies of linearized axisymmetric adiabatic 
perturbations of a purely rotating flow. The criterion for the 
stability of such a flow is that the quadratic form in £ : 


I = 


"a! 

S atu VS o VS \ 


blnq 


c 2 + 


l Iff' 


d Z 


Vlng 
z 


Mna 

& z btff 


<>•» 


is positive definite everywhere. Here g^ and g z are the tor and z 

components of the apparent gravitational acceleration g ’ (which includes 

2 ~ 

the centrifugal acceleration) C is the circulation vs fl and q is the 
potential density [i.a. the density the fluid would have if brought 

to a standard pressure; thus the difference dq in potential density 
between two neighbouring points is dp - (^) s dp » ~ ( tint )p ds 


if the fluid is chemically homogeneous] . 
definite if the conditions 


A 


1 


1 bC* + Vlnq 
0 3 btff btff 


The form I(£)is positive 


> 0 


and 

A 2 = (Vq , VC 2 ) . (g\ «) > 0 


(3.4) 


(3.5) 


are simultaneously satisfied. In deriving the second condition the 
curl of the momentum equation governing the steady state in’ the form 


g\* Vlnq 


1 be 2 

^ Z 


(3.6) 


was used, where is a unit toroidal vector. These conditions were 
obtained earlier from an approximate analysis by Solberg [48] who 
considered the motion of a displaced parcel of fluid in pressure 
equilibrium with its surroundings. 

The analysis of Holland provides sufficient but not necessary 
conditions for stability. A more recent derivation, which amounts almost 
to the same thing, is presented by Fricke and Smith [49] using a variational 
principle derived by Lynden-Bell and Ostriker [50] . Using the variational 
principle is potentially more powerful, for if it can be shown that the 
eigenfunctions of the linear stability problem form a complete set it 
is a simple matter to show without solving the equations that conditions 
(3.4) and (3.5) are also necessary. 
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In certain special cases conditions (3.4) and (3.5) reduce to 
previously known criteria. When ft » 0, * 0 and equation (3.6) 

implies that Vq and g are parallel. Condition (3.4) then reduces to 

the well known requirement that for stability the potential density must 
decrease upwards; for a compressible atmosphere this is the condition 
for convective stability [51, 52]. For rotation independent of z 
in the absence of gravity, A 2 is again zero; if the flow is incompressible 
q ■ p and (3.4) is just the condition (3.2). 

In the more general barociinic case is not zero. As pointed 

out by Hsiland [47] , since Vq and glare not parallel there is a region 

of directions for a displacement € within which (€»g v ) (€,?<!) < 0, so 

that such a displacement would liberate energy via the buoyancy. 2 
Similarly there is a region of directions for £ satisfying (C.trr) (£,VC ; < 0 

within which rotational energy could be liberated. The inequalities 
(3.4) and (3.5) are simply the conditions that these two regions do 
not intersect, and that differencial rotation and the combined effect of 
gravity and the centrifugal force on the density stratification always act 
in opposition. They are satisfied if the circulation increases away from 
the rotation axis in surfaces of constant q [53], or equivalently if q 
increases in the direction of apparent gravity in surfaces of constant C • 

In the discussion above perturbations <5$ in the gravitational 
potential have been ignored. This is a good approximation for disturbances 
of short wavelength, but is not necessarily so otherwise. Taking <$5 
into consideration does not alter the criterion for convective stability 
in the absence of rotation because here the neutral modes of the marginal 
state generate no density perturbations and so 6$ * 0 [54,55]. But 
in a rotating configuration one would not expect that to be the case 
if q is not constant in the equilibrium configuration. Fricke and 

Smith [49] have shown that the perturbation 6$ is never stabilizing, 
so criteria (3.4) and (3.5) are not even sufficient for stability to 
large scale axisymmetric disturbances, and Fricke [56] has constructed 
a cylindrically symmetric example in which perturbations in gravitational 
potential are destabilizing. 

Very different stability criteria are obtained when the perturbations 
are no longer constrained to be adiabatic. Now thermal diffusion can 
weaken the buoyancy forces responsible for stabilizing adiabatic 
disturbances, and provided there is some driving force, sufficiently slow 
motion cannot be prevented. This was pointed out by Yih [57] who 
demonstrated that in the absence of gravity the steady flow of a low 
Prandtl number fluid between two rotating coaxial cylinders can be 
unstable if the Rayleigh criterion (3.1) is not satisfied, even though 
(3.2) is. He showed also that under some circumstances viscosity can 

destabilize an otherwise stable flow if 4r < 0. 

am 

A consistent analysis for barociinic flows has not been 
published, though one can guess at the stability criteria. If (3^4) 
and (3.5) are not simultaneously satisfied there are directions in 
which both buoyancy and rotation tend to drive a displaced fluid element. 
Thermal diffusion acts to reduce buoyancy, but does not change its 
sign, so instability is likely to occur. Again, perturbations in the 
gravitational potential are being ignored. If (3.4) and (3.5) are 
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satisfied, then the criterion for stability can be estimated by setting 
the buoyancy terms in (3.3) to zero [cf 53] : 1(6) is then positive 

definite if and only if 

W- > 0 “* • »• »•’> 

These conditions were first obtained by Goldreich and Schubert [58] 
and Fricke [56] , using a plane wave approximation to the eigenfunctions 
of the linear stability analysis. It was suggested that the conditions 
were necessary and sufficient for stability. In addition to Rayleigh’s 
criterion (3.1) it is also required that ft be independent of z if the 
star is to be stable. This condition is automatically satisfied if 
Vq * 0 in the equilibrium configuration, but not otherwise : Vq was 

set to zero in (3.3) because it was argued that thermal diffusion destroys 
buoyancy in a sufficiently slowly moving flow, even though it is not 
zero in the equilibrium state. 

The growth of the instabilities discussed by Yih, Goldreich and 
Schubert, and Fricke is controlled by thermal diffusion, though the energy 
source is the rotation. The dynamics is essentially the same as that 

of the Eddington-Sweet circulation. Nevertheless the growth time is 

much shorter than the Eddington-Sweet circulation time because these 
modes can occur with very short wavelength; the reduction of the growth 
rates by buoyancy cannot be eliminated entirely, however, because at 
extremely short wavelengths viscous forces come into play. It appears, 
therefore, that criteria (3.7) are likely to apply even to steady 
configurations with meridional circulation. 

The instability of Goldreich and Schubert and Fricke which seems 

to occur when has a profound implication concerning rotating stars 

since it appears that there are no steady solutions to the equations 
governing the structure of a nondegenerate star for which ft is independent 
of z, any steady state that might exist must be unstable. Although 
no proof of the instability has yet been found, it seems likely that „ . _ 

the result is correct. The arguments supporting it predict also the 
instability demonstrated by Yih, which is hardly different, and can also 
be used to predict similar diffusively controlled instabilities such as salt 
fingers [59] , whose existence has been well established [60] . 

The discussion so far has been concerned solely with axisymmetric 
instabilities. Less progress has been made with nonaxisymmetric 
instabilities because they are more difficult to analyse. Rayleigh [61] - 

showed that steady shear flow of a homogeneous incompressible inviscid 
fluid in which the motion is either rectilinear or pure rotation is stable 
to two dimensional infinitesimal perturbations in the plane of the shear 
unless the vorticity somewhere has a turning point. Fj0rtoft [62] showed 
that this turning point must be a maximum in the magnitude of the vorticity. 
General necessary and sufficient conditions have not been obtained [63] 
and it appears that each flow must be analysed separately. However, 
the results suggest that a concentration of vorticity that isn’t bounded 
too closely by rigid walls is unstable. 

Vertical density stratification can stabilize a horizontal flow 
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with velocity u(c) where t is a vertical coordinate, provided there is no 
thermal diffusion and the locally defined Richardson number 

Ri « - g'f| / (§| 2 . 0.8) 

is everywhere greater than ^ [64,653 * Thermal diffusion reduces the 

influence of a gravitationally stable density distribution, just as in 
Yih* s problem, and instabilities can develop on a diffusion timescale [66] . 
Guesses at how these results generalize to more complicated flows must 
be made with caution because the dynamics of shear instability is not 
thoroughly understood. One is warned, for example, by Huppert's [67] 
result that a gravitationally stable density gradient can destabilize 
an otherwise stable shear. 

It is argued by Zahn [68] that shear instabilities provide an 
important angular momentum mixing mechanism in rotating stars. Zahn 
pointed out that all flows at very high Reynolds number appear to be 
unstable, and shear in isentropic surfaces in particular will not be 
stabilized by a gravitationally stable density gradient. The important 
distinction between flows characterized by a high but finite Reynolds 
number and inviscid flows is that the former can generate vorticity 
with viscous stresses, though in stars this must occur in thin boundary 
layers if a resulting instability is to be generated in a time short 
enough to be interesting. In terrestrial conditions it is usually 
against a rigid wall that vorticity is so created. This cannot be the 
case in stars, but internal boundary layers such as those discussed by 
Mestal [33] in the vicinity of a near discontinuity in chemical composition 
can arise, and coupling with other motions such as convection may be 
important. If shear instabilities do commonly occur, then, as Zahn 
stresses, angular momentum transport resulting from them will generally 
dominate any transport by the diffusivily controlled instabilities 
discussed above. 

EVOLUTION OF THE ROTATION OF STARS 

Since it appears that no stable steady state of a nondegenerate 
star exists, stars must evolve to states of time dependent motion so 
long as there is nuclear energy enough to keep them on the main sequence. 
The state to which a star evolves finally, however, after all of its 
available nuclear fuel has been exhausted must be a state of minimum 
energy, which is rigid body rotation. This assumes, of course, that 
an equilibrium state is available to it and that it has not collapsed 
into a black hole. 

The angular momentum distribution in a main sequence star is 
therefore determined by the competing transports due to large scale 
meridional circulation and the smaller scale probably turbulent motions 
arising from the rotational instabilities discussed in the previous 
section. In addition, angular momentum is transported by convection 
and by large scale oscillations, if they are present. 

To what distribution of angular momentum does a star evolve ? 

The only star for which we have direct evidence is the sun. Observations 
of the surface show that the equator rotates faster than the poles by 
about 20%, the precise amount varying somewhat with time, and that the 
rotation of features directly related to variations in the magnetic field 
is greater than the results from Doppler measurements. Magnetic field 
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variations are produced by fluid motions beneath the photosphere, where 
the kinetic energy density exceeds the magnetic energy density, so the 
latter observation suggests that the angular velocity of the sun at least 
near the surface increases with depth (though Yoshimura [69] has argued 
against this interpretation). This is not inconsistent with Dicke's [12] 
conclusion that his measured oblateness of the radiant intensity of the 
solar disk results from an oblateness in the gravitational potential 
caused by rapid rotation of the interior. The situation has been reviewed 
recently by Gilman [70] . 

Little is known about angular momentum transport by motions 
resulting from the instabilities to which a star is subject. It is 
usually assumed that turbulence in a convection zone, for example, acts 
on large scale motion in the manner of molecular viscosity, forcing 
the zone towards a state of rigid rotation. Durney [71,72] , Gilman 
[733 and Gierasch [74], for example, have studied models of the solar 
convection zone in which meridional circulation is acted upon by such 
a turbulent viscosity with a view to explaining the observed differential 
rotation. The circulation is driven by the horizontal pressure gradients 
resulting from rotational distortion, just as in radiative zones. In 
Durney T s models an assumed preferred rotational inhibition of the convective 
heat flux in polar regions contributed further to these gradients [cf 753. 

The models can produce equatorial acceleration at the surface, though 

the dependence of angular velocity on depth is not in all cases that 

suggested by the observations, and in all cases but Durney T s latest 

model the variations in surface energy flux from pole to equator are 

too great to be consistent with observation. The models are not consistently 

coupled to the radiative interior and, as Gierasch has emphasized, this 

may have led to serious errors in the results. 

Another possible source of error is the assumption that the 
Reynolds stresses act like a scalar viscosity. This could be a reasonable 
approximation only if the turbulence were isotropic, but since buoyancy 
and rotation both impose preferred directions on the forces driving the 
flow, anisotropy in the Reynolds stresses is inevitable. Wasiutyhski^63 
and Biermann [77] pointed out that anisotropy arising from buoyancy would 
drive nonuniform rotation in stars, and Kippenhahn [78], Cocke [79] 
and KBhler [80] subsequently employed simple anisotropic stress tensors 
to compute differentially rotating solar envelopes. Models with equators 
rotating faster than the poles can be constructed, with an appropriate 
choice of turbulent stresses, but in the latest ones [80] the angular 
velocity decreases with depth. The anisotropy induced by the rotation 
itself has been discussed by Gough and Lynden-Bell [813 who argued that 
there is a weak tendency for turbulence to force motion on a larger scale 
towards a state of no vorticity. Experiments were performed to support 
their arguments, but later Str ittmatter , Illingworth and Freeman [823 
showed that they had been misinterpreted; the vorticity expulsion hypothesis 
remains unconfirmed. If it is true, however, and if all other transport 
mechanisms are ignored, the Reynolds stresses would act to distribute 
the angular velocity in the solar convection zone such that it decreases 
with latitude and increases with depth. The ideas that led to vorticity 
expulsion liave not been developed enough for Reynolds stresses to be 
calculated. 

Busse [83,84] has studied the dynamics of the convection itself 
by expanding the solutions of the equations of motion about the marginally 
stable state. The conditions are highly idealized and so the results 



25 


are not directly applicable to solar convection, though it is hoped 
that some aspects of the solar dynamics are revealed* In the models, 
the motion develops an anisotropy which leads to increasing surface angular 
velocity from pole to equator* 

A completely consistent picture of the rotation of the solar 
convection zone is not yet at hand, but considerable progress has 
clearly been made in the last decade* 

In computations of the structure of rotating stars, convective 
cores are invariably assumed to be constrained by large isotropic 
turbulent stresses. Tayler [43] has discussed the effect of the anisotropy 
that must be present and concluded that for practical purposes it is 
probably not necessary to take into account the anisotropy in energy 
flow due to convection, but that one must, of course, worry about the 
distribution of angular momentum. 

Angular momentum transport resulting from the development of 
rotationally driven instabilities has been discussed less. Kippenhahn 
[85] pointed out that the transport of angular momentum by the diffusive 
modes discovered by Yih, Goldreich, Schubert and Fricke would lead to 
the formation of a shear which would become unstable to nonaxisymoetric 
instabilities* He argued that this would limit the development of the 
diffusive modes, and by a simple-minded argument estimated the time 
scale for gross readjustment of a star's angular momentum to be at 
least the thermal diffusion (Kelvin-Helmholtz) time t^. James and 

Kahn [86,87] looked at the development of the instability more carefully 
and concluded that the shear would be concentrated in narrow corridors. 

Thus there would be separate regions of diffusively controlled motion 
and dynamically unstable shear. This is not dissimilar to the 
separate layers of fingers and convection in thermohaline convection 
[60,88,89]. James and Kahn estimated that the time for global readjustment 
of angular momentum is as great as the Eddington-Sweet circulation time 
Tgg. It appears that the timescale obtained was as long as t E g because 

the possibility that motion in the diffusive region be on the small 
viscously controlled scale realized by salt fingers was overlooked* 

If' fingering does occur, angular momentum should be transported globally 
in a time shorter than t^g and possibly However such argument 

by analogy may be too naive* The vectorial nature of angular momentum 
and the complicated geometry may act such that this comparatively simple 
layered structure cannot persist, and the growth time of the large scale 
diffusive instabilities, namely t^g, may then be the shortest time in 

which global angular momentum redistribution can occur. 

Throughout this discussion I have had in mind that the angular 
momentum of a star is conserved* This is not always so. For example, 
Kraft [51 has shown that the rotation rates of solar- type main sequence 
stars decrease with age. This appears to be a result of angular momentum 
loss at the surface associated with mass loss in stellar winds £ 90 }* 
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Stars like the sun have surrounding their photospheres high temperature 
coronae which are thought to be maintained by the dissipation of energy 
from mechanical waves which originated in their convection zones below* 

The coronal pressure is too great merely to hold up the outer atmosphere 
hydrostatically [913 , and gas streams outwards exerting a torque on the 
subphotospheric regions of the star via a magnetic field* This process 
has been studied by Ferraro and Bhatia [92], Modisette [93], Weber and 
Davis [94] and Mestel [95, 96], and for the sun leads to significant deceleration 
in a time comparable with its age* The mass lost is a negligible fraction 
of the stellar mass. 

The depth of the convection zones of main sequence stellar models 
decreases abruptly with mass at about M * lj M e ; and it has been 
suggested [e.g. 90] that only deep convection zones can support strong 
coronae and stellar winds. Consequently it is plausible that only stars 
at the low mass end of the main sequence have been significantly decelerated 
on the main sequence, which would explain at least in part the strange 
dependence of angular velocity on mass mentioned in the introduction.. 

However Kraft 1 s [53 observations of solar-type stars in the Pleiades 
group, which have arrived on the main sequence only recently, also show 
an abrupt change in surface rotation rate* It is possible that this 
is a result of very different angular momentum distributions inside upper 
and lower main sequence stars. Alternatively it reflects different 
total angular momentum, which must either have been present at the time 
of formation of the stars or is a product of Kelvin-Helmholtz contraction 
to the main sequence. During part of the gravitational collapse phase 
all stars have extensive convective envelopes, and many are probably 
convective throughout, so conditions seem ripe for something like the 
stellar wind deceleration mechanism to operate very efficiently. There 
are no reliable estimates of the angular momentum that should be lost 
by this mechanism, but since the contraction time decreases with stellar 
mass it is plausible that the low mass stars have been slowed down the 
most. It is perhaps worth mentioning at this point that vorticity 
expulsion in a fully convective star implies angular momentum expulsion 
[81] , and estimates of the expulsion rate during gravitational contraction 
based on analogy with magnetic field expulsion [973 leads to a sharp 
break in the main sequence rotation rate at just about the place where 
it is observed. 

The continual extraction of angular momentum from the sun adds 
yet further complexity to the problem of determining its present state 
of rotation. The magnetic field in the solar atmosphere is too weak 
to be dynamically important much below the photosphere, and so can 
transmit a torque only to the outer regions of the sun. Dicke [12] 
claimed that therefore only the outer convection zone has been decelerated, 
because viscous diffusion would take about 1013 yr to affect the core, 

which is much longer than the age of the sun (about 5 x 10 yr) . The 
argument is invalid, for aside from the existence of the angular momentum 
transport processes more efficient than viscosity, discussed above, the 
changing centrifugal force field near the surface induces pressure 
inbalance which drives currents like the Eddington-Sweet circulation, 
causing the entire sun to respond faster than it would by viscous diffusion 
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alone* this process is called spin-down, and has been reviewed recently 
by Benton and Clark £981 •* ** occurs in homogeneous incompressible 

fluids in rigid containers, angular momentum being exchanged between 
fluid and container in thin viscous boundary layers which were first 
discussed by Ekman [99] ; Einstein [100] , Prandt [101] and Bondi and 
Lyttleton [102] studied the problem further* The sun is not in a rigid 
container, so the details must be somewhat different, but estimates of 
the time to. spin down the solar core based on this analysis and postulating 
a turbulent Ekman layer at the base of thq convection zone [103] , or 
on an analysis treating the convection zone as a rigid porous medium 
[104] , suggest values considerably less than the solar age. However, 
aside from the details of the coupling of the convection zone to the 
radiative interior, these estimates must be treated cautiously because 
they were obtained without taking the stratification of the sun into 
account • 


For a while there was some dispute over the effect of gravitationally 
stable thermal stratification [98] but now it appears, as one would 
naively expect, that thermal diffusion can in general remove any buoyancy 
tending to prevent the spin-down currents* The timescale of the resulting 
flow depends on the properties of the fluid, and linearized theory suggests 
that it is certainly bounded above by [98,105,106] • Analogous results 

hold when stratification is by chemical composition ,but in stars the 
spin-down timescales are then much greater because particle diffusion 
acts more slowly than thermal diffusion. Experiments are in general 
agreement with the theoretical findings [107]* Compressibility appears 
to make no qualitative difference to the conclusions [108] 

The stability and spin-down arguments all lead one to the conclusion 
that global angular momentum readjustments occur in stars in a time 
which is at most the Eddington-Sweet time, but is probably much less* 

In any case, if the sun arrived on the main sequence rotating ten or 
twenty times faster than its surface is now, as Dicke [12,18] believes, 
angular momentum will already have been redistributed at least throughout 
its chemically homogeneous envelope* Approximate calculations by 
Sakurai [38] considering the evolution of the large scale circulation 
in an initially uniformly rotating solar model are consistent with this 
conclusion* If, on the other hand, the entire sun has always rotated 
on the main sequence at a rate comparable with its present surface value 
the issue is less dear* It can be argued that gravitationally stable 
gradients in chemical composition generated in the core by nuclear 
reactions might have prevented the core from spinning down* This is 
possible, though it has been suggested that the core is periodically 
mixed [109] , which would render this conclusion less likely. 

From a fluid dynamical point of view it is therefore quite 
uncertain what the angular momentum history of the sun has been, though 
in balance the evidence favours that a global readjustment has taken 
place on the main sequence* This does not necessarily imply that the 
solar core is not in a state of comparatively rapid rotation, however, 
for advection of. angular momentum by meridional currents generally leads 


The problem usually studied is one where the external torque increases 
the rotation rate. This is called spin-up* 
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to higher angular velocity near the rotation axis [cf 38] . This is 
opposed by turbulent mixing which may tend to make the angular velocity 
more uniform. 

Observational evidence that may have some bearing on the degree 
of mixing that has taken place concerns the abundances of lithium and 
beryllium* Beryllium is destroyed by nuclear reactions under conditions 
prevailing within the inner 50% by radius of the sun* The agreement 
of the solar surface abundance of beryllium with the relative abundance 
in the chondritic meteorites [e.g. 110] might lead one to believe that 
significant mixing of the sun f s surface layers with the interior has 
not occurred. On the other hand lithium, which is destroyed inside 
a radius of about 0.6 R® , is almost absent in the surface* By the 
same argument it appears that mixing down to this level has occurred, 
and since it is below the bottom of the convection zone (which is about 
0.8 R 0 ) rotationally induced mixing of some kind has presumably taken 
place. There is no clear fluid dynamical reason why partial mixing 
to this degree should occur . Dicke [111] has concluded that ‘this 
apparent lack of complete mixing implies that the solar core is not coupled 
to the surface, and is therefore still rotating at its initial rapid 
rate. Contrarily, Goldreich and Schubert [58] argued that incomplete 
mixing is indicative of slow rotation throughout the entire main sequence 
history of the sun, and that the core is therefore unlikely to be rotating 
rapidly. The problem is unresolved, and is complicated by the possibility 
that element separation within the fluid obscures the implications surface 
abundance measurements have on the. degree of mixing that has occurred [112] 

Finally, let us return to the solar oblateness measurements. 

Dicke and Goldenberg [16,17] have measured the oblateness of only the 
radiative intensity distribution in the solar image, and interpreted 
this as a direct consequence of a similar oblateness in the gravitational 
equipotentials . This interpretation has been questioned [113-116] 
and it was suggested in particular that effects of radiative transfer 
and intensity inhomogeneities in the diffuse solar atmosphere may have 
been responsible for the apparent oblateness. Recently Hill and Stebbins 
[117] have measured the shape of the sun, paying careful attention to 
surface nonuniformities [118] . They found time varying brightness 

fluctuations sufficient to account for Dicke and Goldenberg' s 
measurement without recourse to assuming any deviation from the spherical 
shape, and during a period when the brightness was nearly uniform they 
inferred an oblateness no greater than what one would expect if the sun 
were rotating at approximately the surface angular velocity throughout* 

No doubt the controversy is not yet over, but at present it cannot be 
said that the sun provides convincing evidence against the theory of 
general relativity. 
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